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Abstract 



This paper presents some asymptotic results for statistics of Brownian semi-station- 
ary (BSS) processes. More precisely, we consider power variations of BSS processes, 
which are based on high frequency (possibly higher order) differences of the BSS 
model. We review the limit theory discussed in [4, 5] and present some new connections 
to fractional diffusion models. We apply our probabilistic results to construct a family 
of estimators for the smoothness parameter of the BSS process. In this context we 
develop estimates with gaps, which allow to obtain a valid central limit theorem for 
the critical region. Finally, we apply our statistical theory to turbulence data. 
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1 Introduction 

In this paper we study the probabilistic limit behaviour of (realised) power variation, based 
on higher order differences, in relation to the class of Brownian semi- stationary (BSS) 
processes. This class, which was introduced in [12], consists of the processes (Y t )t<=m. that 
are defined by 
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where /j is a constant, W is a Brownian measure on M., g and q are nonnegative deter- 
ministic weight functions on M, with g(t) = q(t) = for t < 0, and a and a are cadlag 
processes. When (a, a) is stationary and independent of W, then Y is stationary, which 
motivates the name Brownian semi-stationary process. 

In the context of stochastic modelling in turbulence the process a embodies the inter- 
mittency of the dynamics. For detailed discussion of BSS and the more general concept 
of tempo-spatial ambit processes see [7, 8, 9, 10, 11, 12]. Such processes are, in particu- 
lar, able to reproduce the key stylized features of turbulence data, such as homogeneity, 
stationarity, skewness, and isotropy. In general, BSS processes are not semimartingales 
(see the discussion in Section 2.2.1 for more details). In consequence, various important 
asymptotic techniques developed for semimartingales, see for instance [6, 19, 28], such 
as the calculation of quadratic variation by Ito calculus algebra and those of multipower 
variation, do not apply or suffice in the BSS setting. 

This study includes a review and some extensions of the theory developed in [4, 5]. 
In this paper we will mainly consider BSS processes without drift (i.e. a = 0), which we 
denote by X. First, let us recall that a properly normalized version of power variation 

[t/A n ] 

22 AfX := X iAn - X (i _ 1)An , 

i=i 

constitute a consistent estimator of the stochastic quantity JjJ \a s \ p ds. Such quantities 
play a similar role in turbulence as in finance: they represent a variability measure of 
the process X. The required normalization depends on the parameter a, which describes 
the behaviour of the weight function g near 0. More precisely, we consider functions 
g with g{x) ~ x a as x J, 0. We call a the smoothness parameter as it describes the 
smoothness properties of X and it corresponds to Kolmogorov's scaling law in turbulence. 
The parameter a also controls the weak limit regime and the rate of convergence for the 
standardized version of power variation. For a subrange of parameters a power variations 

— 1/2 

turn out to be asymptotically mixed normal with convergence rate A n , while other a's 
lead to Rosenblatt type of limits and slower convergence rates. We refer to [31] for the 
corresponding results in the Gaussian setting. 

The main objective of the paper is to use the asymptotic theory for power variations 
to construct efficient estimators of the smoothness parameter appearing in the weight 
function g. It turns out that it is preferable to consider higher order increments when 
computing power variation statistics. In our setting using higher order differences has 
the following crucial advantages: (a) The asymptotic mixed normality of power variations 
becomes valid for all smoothness parameters of interest (a known effect for Gaussian 
models; see e.g. [18]), (b) The power variations are more robust to the presence of smooth 
drift functions. We will review these properties in Section 3 and present some new results. 
In the next step, we will apply the asymptotic theory for power variations to construct 
a consistent estimator of the parameter a. Our main tools are realised variation ratios, 
which compare power variations at two different sampling frequencies (A n and 2A n ). In 
the critical regime we will use estimators with gaps to obtain rate optimal estimators of 
a. This type of estimators were studied by Lang and Roueff [22] in the purely Gaussian 
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framework. 

Finally, we will apply our statistical methods to a turbulence data set obtained from 
one-point measurements of the longitudinal component of the wind velocity in the atmo- 
spheric boundary layer. Our estimator lies around a = —1/6, a parameter value that 
corresponds to the celebrated Kolmogorov's 5/3-law. Furthermore, our statistics show a 
rather stable behaviour when the power p is ranging between 1.5 and 2.5. These facts 
demonstrate that BSS processes constitute adequate models for turbulent flows. 

The paper is organized as follows. In Section 2 we introduce common notation, def- 
initions, main assumptions and some probabilistic properties of the procesess under con- 
sideration; in particular, we establish new connections between the BSS processes and 
certain fractional diffusions processes. In Section 3 we develop the limit results, in Section 
4 we construct estimators for the smoothness parameter of the weight function g, study 
their properties and in Section 5 we apply them to turbulence data. 



2 Definitions and first probabilistic properties 

Let (fi, J 7 , (J r t)teR, IP) be a given filtered probability space. We shall consider a BSS 
process (XtjteM., without drift, defined by 

X t = I g{t-s)a s W(ds), t € R, (2.1) 

J —oo 

where W is an (.F^gR-adapted white noise on M, g : M + — > M is a deterministic weight 
function satisfying g £ L 2 (M + ), and a is an (.F^gR-adapted cadlag intermittency process. 
By an (J r t)-adapted white noise we understand a zero-mean Gaussian random measure on 
Borel sets A C R such that m{A) < oo, with covariance 

E[W{A)W(B)] = m(A n B), 

where m is the Lebesgue measure. Moreover, if A C [t, oo) then W (A) is independent of J~t 
and if A C (— oo,i] then IF(A) is JVmeasurable. Such a random measure is a particular 
case of a martingale measure. Note that for a € M the process (W([a, t])) t>a is a standard 
Brownian motion. Integrals with respect to martingale measures are an extension of Ito 
integrals, see [32] for their definition. 

In order to guarantee the a.s. finiteness of the integral, we furthermore assume that 

/ g 2 (t - s)a 2 s ds < oo a.s., (2.2) 

for any t & R. This condition is obviously satisfied when (cx^tgR is a stationary process 
with a finite second moment (recall that g E L 2 (IR + )). We assume that the underlying 
observations are 

X iAn , i = 0,...,[t/A n ], 

with t > being fixed and A n — > 0. This type of sampling is customarily called infill 
asymptotics or high frequency data. Our class of statistics is based upon higher order 
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differences of the BSS process X. We briefly recall the definition: For any k € N, the fc-th 
order difference A^X at frequency vA n , where d£N, and at stage % > vk is defined by 



3=0 

However, when v = 1 we usually write A™ fc X instead of A™£x. For example, 

A"i^ = X iAn - X(j_ 1)Aji and A™ 2 X = X iAn - 2X ( j_ 1)An + X (i _ 2 )A„- 

In this paper we will restrict our attention to the higher order differences with k > 1, 
although other filters may be considered. In order to define power variations of the process 
X we need to introduce some further notation. We consider a centered stationary Gaussian 
process G = (GtjteR, which we will call the Gaussian core of X, that is given as 



G 



:= I g(t- s)W(ds), tel. (2.3) 

J — OO 



Note that Gt < oo since g G L The correlation kernel r of G is given via 

fn°° q(u)q(u + t)du 
r(t) = Jo ^ l K t > 0. 

\\9\\l 2 (k+) 

Most crucial quantity for the asymptotic theory is the variogram R, i.e. 

R(t) := E[(G t+s - G s ) 2 } = 2|| 5 ||| 2(R+) (1 - r(t)). (2.4) 

We introduce two closely related power variations based on increments A"^ of order k at 
frequency vA n as 

[t/A n ] 

V(X,p,k,v;A n ) t := ^ \A?£X\*>, (2.5) 

i=vk 

V(X,p,k,v;A n ) t := A n T k (vA n )- p V(X,p,k,v;A n ) t , (2.6) 



where p is a positive power and Tk(vA n ) := yE[| A"^G| 2 ]. To determine the limiting 

behaviour of the statistics V(X, p, k, v; A n )t, we need to introduce a set of assumptions, 
which we discuss in the next subsection. 



2.1 Main assumptions 

We start with various conditions on the weight function g : M + — > R. Below, all functions 
Lf : R + — > R, indexed by a given mapping /, are continuous and slowly varying at 0, i.e. 
lim x ^o Lf(tx)/Lf(x) = 1 for any t > 0. Furthermore, the function /( m ) denotes the m-th 
derivative of / and a denotes a number in (— |, 0) U (0, \). 
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(Al): It holds that 

(i) g(x) = x a L g (x). 

(ii) g( k \x) = x a ~ k L g ( k ){x) and, for any e > 0, we have G L 2 ((e,oo)). Further- 
more, \g^\ is non-increasing on the interval (a, oo) for some a > 0. 



(iii) For any t > 



oo 



g^{s)\ 2 a 2 _ s ds <oo. (2.7) 



We remark that in the case of a = the assumptions (Al)(i) and (Al)(ii) are never satis- 
fied simultaneously, hence we excluded this case. The next set of assumptions deals with 
the variogram R. 

(A2): For the smoothness parameter a from (Al) it holds that 

(i) R(x) = x 2a+l L R (x). 

(ii) R( 2k \x) = x 2a ' 2k+1 L R(2k) {x). 

(iii) There exists a b G (0, 1) such that 

l r(^) (y) 



lim sup sup 

xiO y£[x,x b ] 



Lr{x) 



< OO. 



This set of assumptions is standard in the literature; see e.g. [5, 17]. These conditions are 
required to develop the asymptotic theory for power variation of higher order differences of 
the Gaussian process G defined at (2.3). Although the variogram R is uniquely determined 
by the weight function g, assumption (Al) does not imply (A2) in general. However, when 
the slowly varying function L g is smooth enough and satisfies lim^Q L g {x) = c S (0,oo), 
the condition (Al)(i) would naturally imply all other conditions from (Al) and (A2) ex- 
cept (Al) (iii) . Finally, we present an assumption on the smoothness of the process a, 
which is required for the proof of the central limit theorems. 

(A3-7): For any q > 0, it holds that 

E[\a t -a s \i]<C q \t-s\^ (2.8) 

for some 7 > and C q > 0. In particular, condition (2.8) implies by the Kolmogorov's 
criteria that the process a has /3-H61der continuous paths for all (3 £ (0,7). 



The assumptions (Al) and (A2) imply certain probabilistic properties of the BSS process 
X, which we explain in the next subsection. 
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2.2 Some probabilistic properties 

This subsection is devoted to probabilistic properties of the processes G and X, which are 
direct consequences of assumptions (Al) and (A2). 

2.2.1 Is a BSS process a semimartingale? 

After defining the class of BSS processes one naturally asks if this is a subclass of con- 
tinuous semimartingales. We start by exploring this question for the Gaussian core G 
defined at (2.3), as it directly influences the fine structure of the process X. Observing 
the decomposition 

G t+ A-G t = f +A g(t + A-s)W(ds)+ f {g(t + A — s) — g(t — s)}W(ds), 

Jt J -oo 

we obtain by formal differentiation 

dG t = g{0+)dW(t) + ( I g'{t - s)W{ds))dt, 

^ J —oo ' 

where we use the convention g' = g^ 1 '. Indeed, the Gaussian process G is an ltd semi- 
martingale when <?(0+) < oo and g' E L 2 (M + ) and this property also transfers to the BSS 
process X under mild assumptions. However, this situation is less interesting from the 
theoretical point of view, since the asymptotic behaviour of power variation of continuous 
semimartingales is rather well-understood; we refer to [6, 19, 28] for more details. The 
assumption (Al)(i) implies that 

g€L 2 (R + ) but g' L 2 (M+), 

because the derivative of g(x) = x a L g {x) is not square integrable near for a £ (-1,0) U 
(O,!)- It can be shown, see [13], that the conditions g{0+) < 00 and g' G L 2 (M + ) are 
also necessary conditions for G to be a semimartingale. Hence, the process G, and so the 
process X (unless a = 0), is not a semimartingale. 

2.2.2 Asymptotic correlation structure and some consequences 

Now, we turn our attention to the local behaviour of the Gaussian core G. Assumption 
(A2)(i) suggests that the small scale behaviour of the increments of G is similar to the 
small scale behaviour of the increments of B H , where B is a fractional Brownian motion 
with Hurst parameter H = a + 1/2 € (0, 1). This connection can be formalized as follows: 
Let rfc >n be the correlation structure of higher order increments A™ fc G, i.e. 

r k ,nU) ■= corr(A" fc G, A™ +j k G), j > 0. (2.9) 
Then the polarization argument and condition (A2)(i) imply the convergence 



lim r k ,n{j) = Pk{j), 

n— >oo 



(2.10) 
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where pk is the correlation function of the fc-th order increments A^B that are defined 
by 



i-j k > 2 > 



3=0 

and H = a + 1/2. For instance, for k = 1 it obviously holds by assumption (A2)(i) that 
R((j + l)A n ) - 2R(jA n ) + R((j - l)A n ) 1a 

1\2H ol„-|2H i I- i |2H 



n,„ ) = ■' " - # + ii" - *r + li- l 

where the right side is the correlation function of the fractional Brownian noise. Con- 
sequently, the limit theory (law of large numbers and the central limit theorem) for the 
power variation V(G,p, k; A n ) t is expected to be the same as for V(B H ,p, k; A n ) t with 
H = a + 1/2. The latter is well- understood since the work of [14]. We shortly recall this 
classical result. Let 

f(x) = \x\P-m p , (2.11) 
m p = E[\U\ p ], U~jV(0,1). (2.12) 

Then the function /, which is associated with a centered version of the statistic 

V(B H ,p,k,l;A n ) t , 

exhibits a Hermite expansion of the form 

oo 

f(x) = J2a l H l (x), (2.13) 

1=2 

where (Hi)i>q are Hermite polynomials and a2 7^ 0. The number 2, which is the minimal 
index I with ai 7^ 0, is called the Hermite rank of the function /. Now, the statistic 
V(B H ,p, k; A n )t is asymptotically normal (with a standard A r ^ 2 rate) whenever the 
condition 

00 

^2p k (j) 2 < 00 

3=1 

holds, where the power reflects the Hermite rank of /. A simple computation shows that 
\Pk(j)\ = 0(j 2 ( H ~ h >) as j — > 00. Hence, the preceding condition is satisfied for (a) k = 1 
and H £ (0, 3/4), (b) k > 2 and H £ (0, 1). When k = 1 we obtain the following cases: 

< H < 3/4 : A' 1 ' 2 {V(B H ,p, 1, 1; A n ) t - m p t) A v p W t , 

# = 3/4: (A n log A^ 1 )^ 1 / 2 (V(B H ,p, 1, 1; A n ) t - m p t) A^, 
3/4<tf <1: A 2 n H - 2 (V(B H ,p,l,l;A n ) t -m p t) Ai t , 

where the weak convergence takes place on B([0,T]) equipped with the uniform topology, 
W denotes a Brownian motion and L is a Rosenblatt process (see e.g. [31]). Finally, the 
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constants v p and v p are given by 

oo oo 

^:=^Z!af(l + 2^ Pl (j) i ), 



fl— 1 r OO 

,2 



5 P := 2 lim -J- £ — Pl (j) 2 • ]T Zlaf 

n— ^oo nffn ^ — « ri ^ — « 



logn ^— ' n 
2.2.3 Connection to integral processes 

In the papers [2, 15, 24] the asymptotic theory for power variation of integrals with respect 
to Gaussian processes and related functionals have been developed. A brief comparison of 
the asymptotic results for such integral processes and BSS models shows that the limit 
theory is quite similar. Thus, a natural question appears: Do BSS models exhibit a 
representation as an integral with respect to a Gaussian process? First, we observe that 

X t + C <r s dG s , 



o 



where G is a Gaussian core of X. However, these two processes are indeed related. Let us 
assume for the moment that the process a is deterministic. We introduce the decomposi- 
tion X t = X[ + X'l with 



X[= f g(t - s)a s dW s , 
Jo 

,0 

X' t '= / g(t-s)a s dW s 



Similarly, Gt = G' t + G", where G', G" are defined exactly as X' , X" with a = 1. Applying 
the integration theory developed in [23], we deduce the identity 

J a s dG' s = X' t + jT (f\*8 ~ °u)g'(s - u)W(du)j ds. 

We remark that the integral on the left side is well-defined in the Young sense when the 
assumption (A3-7) is satisfied for some 7 £ (1/2 — a, 1), since the process G' has Holder 
continuous paths of all orders smaller than a + 1/2. A straightforward computation shows 
that the integral on the right side is well-defined in the Ito sense under the same condition 
(A3-7) with 7 e (1/2 - a, 1). 

The second part X" satisfies the differential equation 



dX'l = (^J g'(t - s)a s W(ds)j dt 



and similar identity holds for G" . Notice that the involved Brownian integral is finite due 
to assumption (Al)(iii) applied to k = 1. Summarizing these findings we conclude that 
Xt can be decomposed as 

X t = a s dG s + A t , 
Jo 
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where A is a continuously differentiable process. Thus, the law of large numbers for power 
variations of X and J a s dG s coincide, since the process A is too smooth to affect the limit. 
However, the central limit theorem for the power variation of X may be seriously affected 
by the presence of the drift A. We will study this type of effects in the next section. 



3 Limit theory for power variations 

Before we present the main limit theorems for power variation of the BSS process X, 
let us introduce some further definitions. For elements Y n ,Y from the space of cadlag 
function D([0,T]), we write Y n ^> Y whenever sup[ 0) T] l^t" — 0- We say that a 

sequence of processes Y n converges stably in law to a process Y, where Y is defined on an 
extension (Of, ,F',P') of the original probability (O, T, P), in the space D([0, T]) equipped 
with the uniform topology (Y n Y) if and only if 

lim E[f(Y n )Z] = E'[f(Y)Z] 

n— >oo 

for any bounded and continuous function / : D([0,T]) — > K and any bounded J 7 - meas- 
urable random variable Z. We refer to [1], [20] or [29] for a detailed study of stable 
convergence. Note that stable convergence is a stronger mode of convergence than weak 
convergence, but it is weaker that u.c.p. convergence. 

First, we recall the law of large numbers for the statistic V(X,p, k, v; A n )f. The cor- 
responding result for k = 1 was proved in [4], while the case k = 2 was treated in [5]. The 
extension to a general k > 1 is straightforward, and therefore omitted. 

Theorem 3.1. Assume that the conditions (Al) and (A2) hold. Then we obtain that 

V(X,p,k,v;A n ) t u -^V(X,p) t :=m p f \<j 8 \*da t (3.1) 

J o 

where the power variation V(X, p, k, v; A n )t is defined at (2.6) and the constant m p is 
given by (2.12). 

We remark that the statistic m~ 1 V(X, p, k, v; A n )t is a consistent estimator of the 
stochastic quantity J* * |cr s | p d,s. However, the computation of this statistic requires the 
knowledge of the parameter Tfc(vA n ), which in turn depends on the weight function g. 
Nevertheless, even when g is unknown, Theorem 3.1 may be applied to estimate the 
smoothness parameter a. We will discuss the estimation procedure in the following section. 

Now, we present the central limit theorem associated with the u.c.p. convergence at 
(3.1). More precisely, we state the central limit theorem for the properly standardized 
vector (V(X,p, k, 1; A n ), V(X,p, k, 2; A n )), since it will be required in the next section. 
The proof for k = 1 (resp. k = 2) can be found in [4] (resp. [5]). The extension to the 
general case of k > 1 follows along the lines of the proof in [5]. 
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Theorem 3.2. Assume that the conditions (Al), (A2) hold and (A3-j) is satisfied for 
some 7 G (0, 1) with j(p A 1) > 1/2. If k = 1 we further assume that a £ (— |, 0). Then 
we obtain the stable convergence 

A' 1 / 2 (v(X,p, k, 1; A n ) t - V(X,p) t , V(X,p, k, 2; A n ) t - V(X,p) t ) A J* \a s \ p A p dB s 

(3.2) 

on D 2 ([0, T]) equipped with the uniform topology, where B is a 2-dimensional Brownian 
motion that is defined on an extension of the original probability space (Q, J 7 , P) and is 
independent of F, and the matrix A p = (A/ )i<ij<2 is given by 

Xl 1 = lim A^ 1 vsa(v(B H ,p,k,l;A n ) 1 ), Xf = lim ^Wf^,?, ft, 2; A n )i) 

^ n— >oo \ / r n— >oo \ / 

A* 2 = lim A^covfv^ ,p, fe, 1; A„)i, K^.p, 2; A n )i) . (3.3) 

with B H being a fractional Brownian motion with Hurst parameter H = a + 1/2. 

Remark 3.3. We remark that the constants Ap defined at (3.3) are indeed finite. For 
instance, it holds that 

oo oo 
1=2 j=l 

where (a^)/>2 are Hermite coefficients that appear in the Hermite expansion in (2.13) and 
the correlation function pk(j) is defined at (2.10). The finiteness of Ai, 1 follows directly 
from \pk(j)\ < 1, SfeiPfc(i) 2 < 00 ( see section 2.2.2) and the identity 

oo 

oo > var[|[/H =m 2p -m 2 = ^Uaf, U~tf(0,l). (3.4) 

1=2 

□ 

Remark 3.4. There are various multivariate extensions of Theorem 3.2. We refer to [4, 5] 
for joint stable central limit theorems for the family 

a- 1 ' 2 (^^.M^^-yfi,^),)^, 

where (pj)i<j<d ar e positive powers and {vj)\<j<d are natural numbers. Such joint central 
limit theorems are important, because our estimator of the smoothness parameter a is a 
ratio of power variation statistics compared at two different frequencies (A n and 2A n ). □ 

We remark that taking higher order difference leads to a central limit theorem for 
all values of a in the interval (— 1,0) U (0, ^), while the convergence at (3.2) for k = 1 
holds only when a G (—4,0). This effect of higher order filters is well-known for Gaussian 
processes; see e.g. [18]. Another advantage of higher order differences is its robustness 
to smooth distortions. Let C m ([0,T]) denote the class of functions / that are [m] times 
continuously differentiable with /([ m D being Holder continuous of order m— [m]. We obtain 
the following lemma. 
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Lemma 3.5. Let (-A*)te[o,T] ^ e a stochastic process with paths in C m ([0,T]) for some 
m > 1. Consider the process 

Z t = X t + A t . 

Under the assumptions (Al) and (A2) it holds that 

\V(Z,p, k, 1; A n ) t - V(X,p, k, 1; A n ) t \ = P (|A^-«- 1 /2 jL (A n )r 1 ), (3.5) 

where L : IR + — )• R is a slowly varying function. In particular, the rate on the right hand 
side converges to 0. 

Proof. The condition A G C m ([0,T]) implies that 

\Al k A\ < CA k n Am 

for some constant C > 0. On the other hand we have that 

r k (A n ) = A a n +1 ' 2 L(A n ) 

for some slowly varying function L, due to assumption (A2)(i). Assume that p £ (0,1]. 
Then we deduce 

\V(Z,p, k, 1; A n ) t - V(X,p, k, 1; A n ) t | 

[t/A n ] 

<A n r fc (A n )^ ^ lA^yll^OpdA^T^A,,)- 1 ^), 

which implies the assertion of Lemma 3.5 for p £ (0, 1]. For p > 1 we apply the mean 
value theorem to conclude that 

\V(Z,p, k, 1; A n )i - F(X,p, fc, 1; A n ) t | 

[t/A n ] 

<pA n T k {A n )-v £ (lA^I + IA^XD^IA^I. 

i=k 

Using that 



lV^n\ f . t 

~I ■'0 



[t/A„ 

L c2s 

and | A^ fc A| < CA^ Am with m > 1, we deduce the approximation 

\V(Z,p, k, 1; A n ) t - U(A>, k, 1; A n ) t | = Op(A,^ Am r fc (A n ) _1 ), 
which implies the assertion of Lemma 3.5 for p > 1. □ 

We remark that the degree of robustness to drift processes A S C m ([0, T]) is increasing 
in fc. In fact, for k = [m] + 1 the right side of (3.5) converges to at the fastest rate. Since 
m > 1, the statistics V(Z,p, k, 1; A n )j and V"(X,p, fe, 1; A n )j are asymptotically equivalent. 
They also satisfy the same central limit theorem given (k A m — a — l/2)(p A 1) > 1/2. 
For instance, when A: = 1 and p > 1, the quantity V(Z,p, k, 1; A n )t satisfies Theorem 3.2 
only for a G (-1/2,0). 
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Remark 3.6. As we mentioned in section 2.2.3, when a is deterministic and (A3-7) is 
satisfied for some 7 G (1/2 — a, 1) (this means that a is Holder continuous of order 
7 G (1/2 — a, l)),we have the decomposition 

X t = I (J s dGs + A t =: X t + A t 
Jo 

with A G C 1 ([0, T]). For k = 1 the central limit theorem for the statistics V(X,p, k, 1; A n )t 
holds for a G (—^,0) U (0, |) (see e.g. [2]), which is in line with the discussion in section 
2.2.2. However, when k = 1, Theorem 3.2 holds for the quantity V(X,p, k, 1; A n )t only if 
a G (— 2,0). This follows from the result of Lemma 3.5, although the rate in (3.5) is by 
no means sharp. □ 



4 Estimation of the smoothness parameter 

We apply the asymptotic theory presented in the previous section to construct consistent 
and asymptotically normal estimators of the smoothness parameter a. We consider mainly 
the specification 

g[x) = x a exp(-\x), x > 0, (4.1) 

where A > is another parameter, but some of the results apply also to a general g that 
merely satisfies conditions (Al) and (A2). We use the change-of- frequency statistic based 
on second order increments 

rni5V a x V(X,p,2,2;A n )t , „ 

COF(p, A n ) t = —— , t > 0, 

V{X,p,2,l;A n )t 

with p > 0, as a foundation for our first estimator. To understand the asymptotic behavior 
of COF(p, A n ) t we write 

'T 2 (2A n ) 2 \ p/2 V(X,p,2,2;A n )t 



COF(p,A 



n t 



r 2 (A„) 2 ; F(X,p,2,l;A 



nit 



Under Assumption (A2), we have 

T2(2A n ) 2 = AR{2A n ) - R(AA n ) 2a+1 
r 2 (A„) 2 AR(A n ) - R(2A n ) 

so by Theorem 3.1, the following consistency result is immediate. 

Proposition 4.1. Under the conditions of Theorem 3.1 we have for any p > 0, 

COF(p,A n )*^2^- 

and, consequently, 

a(p, A n ,t) := V(COF(p, A n ) t ) ^> a, (4.2) 

where 

, / \ !og 2 (^) 1 . n 

h p {x) = -, x > 0, 

p 2 

with log 2 standing for the base-2 logarithm. 
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Note that, since (4.2) holds for any p > 0, computing the change-of-frequency statistics 
for different values of p can be used to gauge the robustness of the estimates of a. 

Next we derive a feasible, asymptotically normal test statistic in the case p > 2 (see, 
however, Remark 4.3 below). To this end, let us consider the decomposition 



(2a + l)p s 
2 



A; 1 / 2 (COF(p,A n ) t -2- 

' T 2 (2A n ) 2 \ p/2 f A- 1/2 (V(X,p,2,2; A n ) t - V(X,p) t ) 



r 2 (A n ) 2 ; V V(X,p) t 
V(X,p, 2, 2; A n ) t A- 1/2 (V(X,p, 2, 1; A n ) t - V(X,p) t 
V(X,p,2,l;A n ) t V(X,p) t 

\2\p/2 



(4.3) 



When g is given by (4.1), one can establish the estimate 

A - 1/2 " 22a+1 ) = 0(A " /2_2Q) ' (44) 

which is sharp in the sense the big O cannot be replaced with a little o; we refer to 
section 5.2 in [3] for the justification of (4.4). Thus, due to the inequality \x r — y r \ < 
r(x r ^ 1 V y r ~ 1 )\x — y\, where x, y > and r > 1, the latter term on right-hand side 
(r.h.s.) of (4.3) is asymptotically negligible when a G (— |, 0) U (0, j) and p > 2. Applying 
Theorem 3.2, the former term on r.h.s. of (4.3) converges stably in law to 

(2a+l)p t .1/2 

2 2 (-1,1) J \a s \PA p < dB s 



V(X,p) t 

Invoking the delta method and the basic properties of stable convergence, we arrive at an 
asymptotically normal test statistic given as follows. 

Proposition 4.2. // the conditions of Theorem 3.2 hold and g is given by (4.1) with 
a G (— 0) U (0, \), then for any p > 2, we have 

(HP, A n ,t)-a)V(X,p,2,l;A n ) t ^ ^ 

\h' p (COF(p, A n )i)|COF(p, A n ) t ^m^V(X, 2p, 2, 1; A n )(-1, l)A p (-l, l) T 

Remark 4.3. If p G (^,2), then under the stronger restriction a G (— \,^^) \ {0} it is 
possible to show that the latter term on right-hand side (r.h.s.) of (4.3) is asymptotically 
negligible and, thus, the conclusion of Proposition 4.2 holds. 

Remark 4.4. An alternative method to estimate the smoothness parameter is using the 
realized variation ratio statistic. However, it does not allow the same flexibility over the 
choice of the power p, as the change-of-frequency statistic does. We refer to [5] for more 
details on the realized variation ratio. □ 
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The estimation method described so far has, of course, the drawback that it applies 
only when a € (— |, 0) U (0, j). In the case a 6 [4,5)5 the latter term on the r.h.s. of (4.3) 
is non- negligible due to the sharpness of the estimate (4.4). To cover the critical region 
[I, h), we develop a theory of power variations with gaps. More precisely, we consider for 
v = 1,2 modified power variation statistics of the form 



V(X,p, k,u n ,v; A n ) 



t ■- 



' [t/(«„A„)] 

53 iA^ njfc xr, v = 1, 

i=[k/u n ]+l 

[t/M„)]-i 

V |A n '\ r /91 ,X| P , « = 2, 

1 Ml„ + [tl rl /2],fe I ' ' 

i=[k/u„] + l 

V(X,p,k 

)t '■= UnA n Tk(vA n ) p V(X,p,k,u n ,v; A n )t, 

where u n £ N denotes the size of the gaps. Simply put, the variation V(X,p, k,u n , 1; A„)t 
is computed by taking only every u n -th increment into account, whereas 

V(X,p, k,u n ,2; A n ) t 

is defined similarly, but choosing increments that fall between those that contribute to 
V(X, p, k, u n ; A n )t via the shift [u n /2] in the definition. In what follows, we let u n — > 00 
so that u n A n — > 0. We also assume that n is always large enough so that u n > {Ak + 2)/3 
(to ensure that the definition above makes sense). 

It is straightforward to show that a law of large numbers, a counterpart of Theorem 
3.1, holds also for power variation with gaps. That is, if conditions (Al) and (A2) hold, 
then 

V(X, p, k, u n , v- A n ) t u -^> V(X, p) t . (4.5) 

for v = 1, 2. Next, we state a central limit theorem for power variation with gaps, which is 
actually simpler than Theorem 3.2. This is due to the fact that the widening gaps between 
the increments make them asymptotically uncorrelated. 

Theorem 4.5. Assume that the conditions (Al), (A2) hold and (A3-j) is satisfied for 
some 7 £ (0, 1) with j(p A 1) > 1/2. Moreover, let u n — > 00 so that u n A n — > 0. If k = 1 
we further assume that a G (— J,0). Then we obtain the stable convergence 

(u n A)- 1 ' 2 (V(X,p, k, u n , 1; A n ) t - V{X,p) t , V(X,p, k, u n , 2; A n ) t - V(X,p) t 



St 



m 2p ~ m p J \ a s\ P dB s 



on D 2 ([0, T]) equipped with the uniform topology, where B is a 2- dimensional Brownian 
motion that is defined on an extension of the original probability space (Q, T, P) and is 
independent of T . 



The proof of Theorem 4.5 is very similar to that of Theorem 3.1. We provide a sketch 
of it, primarily to outline the key differences, in appendix A. 
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To construct our second estimator for the smoothness parameter a, we look into the 
properties of the modified change-of-frequency statistic 

nnv( a \ V{X,p,2,u n ,2;A n ) t 
COF(p,A n ,u n ) t - 



V{X,p,2,u n ,l;A n y 



involving power variations with gaps. Due to the law of large numbers (4.5), it is immediate 
that 



h p (COF(p,A n ,u n ) t ) 



u.c.p. 

— > a, 



under conditions (Al) and (A2), where the function h p is defined in Proposition 4.1. The 
key point of using gaps is that, thanks to the slower rate of convergence (u n A n )~ 1 / 2 in 
Theorem 4.5, we may write 

{UnAn)- 1 ' 2 (COF(p, A n , U n ) t - 2^^) 

' r 2 (2A n ) 2 \ p/2 ( K A n )-V2 (V(X, p, 2, u n , 2; A n ) t - V(X, p) t ) 



r 2 (A„) 2 J V V(X,p) t 

_ V(X,p, 2, 2; A n ) t KA n )- 1 / 2 (y(X,p, 2, u n , 1; A n ) t - V(X,p) t ) 
V{X,p,2,u n ,l;A n ) t V(X,p) t 

+ (u n A n )" 



_ 1/2 ( ( T 2 {2A n ? y' 2 , o2a+nP/2 N 



(4.6) 



r 2 (A n )2 

where, provided that the assumptions of Theorem 4.5 are met, the former term on the 
r.h.s. converges stably in law to 



2 i ^ + ym 2p -m2 ^ \a s \PdW s 
V(X,p) t ' 

where W is a one-dimensional standard Brownian motion independent of J 7 . Again, let 
us assume that g is given by (4.1). Using (4.4), we find that the latter term on the r.h.s. 
of (4.6) is asymptotically negligible if and only if the sizes of the gaps satisfy 

n; 1 = o(A^ 1 ). (4.7) 

(When a € [j, \), for instance u n = [A~ K ], where k G (4a — 1, 1) ^ 0, satisfy both u n A n — > 
and (4.7).) Then, we obtain the following result that characterizes an asymptotically 
normal test statistic. 

Proposition 4.6. If the conditions of Theorem 4-5 hold, g is given by (4.1) with a G 
(— |, 0) U (0, o), and (4.7) holds, then for any p > 2, we have 

{h p (COF(p,A n ,u n ) t ) - a)V(X,p,2,l;A n ) t d 



h' p (COF(p, A n ) t )|COF(p, A n ) n /4(1 - m2/m 2p )F(X, 2p, 2, 1; A n ) 
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Figure 1: Gray dots: The spectral density of the measured data estimated using Welsh's 
overlapping segment averages with a Hanning taper, a segment length of 2 million, and an 
overlap of 50 %. Black curve: The spectral density (5.1) with a = —0.165 and A = 0.0884 
estimated using a least squares method in the double logarithmic representation. Only 
data with frequencies below 200 Hz were used. 

Remark 4.7. Obviously, the price of using gaps is that we lose lots of observations. How- 
ever, such a procedure may be optimal if the observations are highly correlated. In fact, 
Lang and Roueff [22] consider the change-of-frequency statistic in the context of Gaussian 
stationary increment processes and show (see [22, Theorems 1-3]) that the estimator with 
gaps achieves the minimax bound, which is given by A 2 "" 1 when a £ [|, |). In the setting 
of BSS processes the estimator COF(p, A n , u n )t attains the minimax bound if we choose 
u n oc A* -4 "; in that case Proposition 4.6 would contain an asymptotic bias. □ 

5 Application to turbulence data 

To illustrate the estimation of the smoothness parameter a, we will consider a data set 
consisting of 20 million one-point measurements of the longitudinal component of the wind 
velocity in the atmospheric boundary layer, 35 m above ground. The measurements were 
performed using a hot-wire anemometer and sampled at 5 kHz with a resolution of 12 bits. 
The time series can be assumed to be essentially stationary, the mean is 8.3 m/s, and the 
standard deviation is 2.3 m/s. We refer to [16] for further details on the data set (the data 
set is called no. 3 therein). The observations are standardised to have zero mean and unit 
variance. 

The spectral density S of a BSS process with g given by (4.1) satisfies 



S(f) = const, x (1 + (2vr//A) 2 )-( 1+a ), 



(5.1) 
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where / denotes the frequency. Thus, for / 3> A, we have S(f) oc f- 2 ( l + a ) . i n the context 
of turbulence, Kolmogorov's 5/3-law [21] states that in the so-called inertial range, the 
spectral density is approximately proportional to /~ 5 / 3 . Putting a = —1/6 reproduces 
the 5/3-law for / 3> A. The weight function g therefore defines an infinitely long inertial 
range bounded from below, in the frequency domain, approximately by A. 

Figure 1 shows the spectral density estimated from the turbulence data. For frequencies 
below approximately 200 Hz, the weight function g provides a good fit with a = —0.165 ~ 
—1/6, in agreement with Kolmogorov's 5/3-law. The inertial range appears to be from 0.1 
Hz to 200 Hz. At higher frequencies (or, equivalently, smaller scales) the model no longer 
describes the spectral density accurately since the dissipation of the fluid's kinetic energy 
into heat causes the spectral density to decay approximately exponentially fast [30]. 

We use the estimator a(p, A n , t) based on the change-of-frequency statistic, given by 
(4.2), to estimate the smoothness parameter a. Figure 2 shows the estimates &{p, A, t) for 
t = 2 000 000,20 000 000 (xl/5000 sec), p = 1,2,4,8, and A = 1,...,2000 (xl/5000 sec). 
Several comments are in order. Unsurprisingly, the dispersion of the estimated values of a 
increases with increasing p. In the case t = 20 000 000, the amount of data available for 
the estimate of a is large enough to ensure a reliable estimate. We will therefore focus on 
the case t = 20 000 000 (right column) . The best agreement with the expected a = — 1/6 is 
found in the case p = 2 when A > 500. This is perhaps not surprising, since p = 2 makes 
a(p,A,t) a second order statistic, and the spectral density, being the Fourier transform 
of the autocovariance function, is itself a second order statistic. Furthermore, A > 500 is 
well within the inertial range where the weight function g has been shown to agree well 
with the data. Beyond the displayed estimates, a(p, A, t) appears fairly robust to changes 
in p when 1.5 < p < 2.5. At small lags, dissipation-scale effects cause the departure from 
a = —1/6 to be expected. The dependency of a(p, A, t) on p for larger p hints at some 
lack of robustness of the estimator. However, this is likely due to insufficiencies of the 
BSS process in describing the data, as simulations (Gaussian with a = 1) display no such 
dependency on p. 
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A Sketch of the proof of Theorem 4.5 

The main tool for the proof of Theorem 4.5 relies on is the blocking technique used in 
[2, 4, 5]. That is, for any £ > 0, we divide the interval [0,t] into blocks of length t . We 
then keep the value of the volatility process a fixed in each of the blocks and establish the 
stable convergence of the power variation. Finally, we pass to the limit I — > oo giving us 
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Figure 2: Gray dots: The estimate a(p,A,t) for t = 2 000 000 (left column) and t = 
20 000 000 (right column) for p = 1, 2, 4, 8 (from top to bottom). The lag A varies from 1 
to 2000. Black line: The value —1/6 corresponding to Kolmogorov's 5/3-law. 
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the limit that appears in the statement of Theorem 4.5. 

More concretely, the blocking technique is based on the decomposition 

(%A); 1/2 (V(X, p, k, u n , 1; A n ) t - V{X,p) t ) = S(£, n)\ + S{n) 2 + S(£, nf t + S(£, n)f, 

where 



[It] 

S(£,n)l :=J2\*U-W\ P ( u n A n)- 1/2 {u n A n T k (A n )- p / 2 £ \^ Un , k G\ p - m p C x 

j=l i£l e (j,n) 
[t/(u n A„)] 

S(n) 2 := M„) 1/2 r t (A„)- p / 2 £ (\^ n , k X\ p - \a (lUn „ k)An A^ k G\ p ), 

i=[k/u n ]+l 
[t/(«„A„)] 

S(£,nf t := (n n A n ) 1/2 r fc (A n )^/ 2 ( £ l^-j^A^Gl* 

i=[fc/u n ]+i 

[ft] 

i =1 ieli(j,n) 

[ft] 

S{t,n)t := m. 



.{UrAnY^^Y^Wij-D/AV ~ [ Ws\ P ds), 

j=l J o 



with Ii(j,n) := [i G N : i > [£;/«„], iu n A n 6 ((j — l)/£,j/£]}. An analogous de- 
composition can be derived also for y(X, p, A;, -u n , 2; A n )i under the same centering and 
normalization. Using methods similar to [5], we may show that for any 5 > 0, 



limlimsupP sup \S(n)f + S(£, n)f + S(£, n)t\ > S = 0. 

>oo n -¥oo M6[0,T] ' 

Thus, only S{£,n)\ contributes to the limit. Let us denote by B 1 and B 2 the components 
of the 2-dimensional Brownian motion B appearing in the statement of Theorem 4.5. By 
Lemma A.l, below, and the properties of stable convergence, we have that 

[ft] 



S(l,n)\ yjm 2 p - ml^2\a {j _i )/e \ p (Bj /e - B\^-_ l y t ), n^oo. (A.l) 

i=i 

When £ —> oo, the r.h.s. of (A.l) converges in probability to 

sj m 2p - m 2 J \a s \ p dB l s . 

Again, analogous statements hold jointly for V(X,p, k, u n , 2; A n )t, but B 1 replaced with 
B 2 therein. 

To complete the argument, it remains to prove the following lemma concerning the 
power variations of the Gaussian core G. 
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Lemma A.l. Assume that the conditions (Al), (A2) hold. Moreover, let u n —> oo so 
that u n A n — > 0. If k = 1 we further assume that a 6 (—^,0). Then we obtain the stable 
convergence 

(«„A); 1/2 (v(G, p,k,u n ,l; A n ) t - m p t, V(G, p, k, u n , 2; A n ) t - m p tj yj m 2p - m 2 B t 
on D 2 ([0,T]) equipped with the uniform topology, where B is as in Theorem 4-5. 

We need to introduce first some concepts that are needed in proof of Lemma A.l. 
Let us denote by H. the closed subspace of L 2 (Q, F, P) spanned by the Gaussian random 
variables Gt, t > 0. Note that T~L is then a Hilbert space consisting of centered Gaussian 
random variables. Thus, the inclusion map 7i L 2 (Q, T, P) can be seen as an isonormal 
Gaussian process, with respect to which we may define multiple Wiener integrals of orders 
1 = 1,2,..., denoted by Ii. Recall that // is a linear isometry T-L &1 — > L 2 (f2, T, P), where 
T-L &1 denotes the /-fold symmetric tensor product of T~L. In particular, I\(Gt) = Gt for all 
t > 0. 

Suppose that U G H is such that \\U\\ H = E[U 2 ] 1/2 = 1 and / is defined by (2.11). 
Then, using the Hermite expansion of / (2.13) and the connection of Hermite polynomials 
and multiple Wiener integrals [25, Proposition 1.1.4], we obtain a chaos expansion 

oo 

f(U) = J2aiHU® 1 ), (A.2) 

1=2 

strongly convergent in L 2 (Q, J 7 , P), where U® 1 stands for the /-fold tensor product of U. 
Convenient sufficient conditions that functionals admitting chaos expansions of the form 
(A.2) converge to a Gaussian law are provided in [2, Theorem 5], building on the results 
of Nualart and Peccati [26], and Peccati and Tudor [27]. We will rely on these conditions 
in the proof, below. One of the conditions involve the so-called contractions of elements 
of U® 1 . To recall the definition, let Vi, V 2 G H® 1 , where I > 2. Then for q = 0, 1, . . . , I - 1, 
the g-th contraction of V\ and V2, which is an element of H 2 ( l ~ q \ is given by 

{V 1 ® q V 2 )(uj' 1 ,...,uj' 2{l _ q) ) 

■= J Vi(u[, . . .,u)i_ q ,u)i, . . . ,uj q )V 2 (uji, . . . ,uj q ,u}[_ q+1 , . . . ,W2 ( i_ (? ))P® 9 (da;i, . . .,dw q ). 

In the case q = the definition reduces to the tensor product V\ ® V%. 
Proof of Lemma A.l. To simplify notation, let us denote for v = 1,2, 
Z v t := (u n A)-^ 2 (V(G,p,k,u n ,v;A n ) t -m p t). 
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Using (A. 2) and the linearity of multiple Wiener integrals, we obtain 

[t/ (u n A n )] /An v\ 

Z} = (u n A n ) l/2 E f(^ Y )+0((u n A n )^) 

r j / iii V \ n ) J 
i=[k/u n \ + l 

oo / [*/(«„ A n )] , x n Y\® 1 \ 

1=2 V i=]k/u„]+l V fcl nj / J 



where the remainder term is uniform in t. Analogously, 

[t/(u n A n )]-l ,aii,Z X\® 1 

£ 

1=2 \ i=[k/u n ]+l 



oo / [t/(u n A n )]-l ,£U,2 X\® 1 \ 

zf = j2iiU(u n A n )^ e ( ) +°(^ A «) 1/2 )- 



To verify the conditions of Theorem 5 of [2], we need to estimate, for any s > t, the 
inner product 

/T?l,l,n j?l,l,n< 

in® 1 

where the remainder term is uniform in I and 



/ F U„ F U»\ = a 2 (t + A\ n + A\ n + 0(u n A n )), 



[t/(u„A n )]-[k/u„]-l 

Aj n = 2u n A n E ([t/(u n A n )} - [k/u n ] - j)r k ,n(unj) 1 , 

3=1 

[t/(unA„)]— [k/un] 
■^l,n = u nA n ^ Cj,n r k,n( u nj) ; 

i=i 

with supj \cj tn \ < [s/(u n A n )] — [t/(u n A n )]. By Assumptions (Al) and (A2), there exists 
(cf. Lemma 1 of [2]) a sequence (r~k{j))j >0 such that \r k ^ n (j)\ < f k {j) for all j > and 
n £ N, with 

oo 

Yr k (j) 2 < oo. 
3=0 

Since / > 2, we have 

[t/(u n A n )]— [k/u n ]— 1 oo oo 

I4nl < c E ^.,nKi) 2 < cE^Ki) 2 < c E f fc0') 2 ' 

j=l j=l j=u n 

oo oo 

I4 n | < (s - t + o(u n A n )) E ^Ki) 2 < c E ^(i) 2 > 

where C > is a constant that depends on s, but not on n nor I. Letting n — > oo we 
obtain now 

?l,l,n r?l,l,n\ 



l\(F^,F^ n ) Hm ^llaft 
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due to the assumption u n — > oo. In the view of (3.4), it is then immediate that 

oo 

lim limsup ZIIIF. 1 ' = 0. 

nwoo n _ >00 ^ 
l=m 

We may establish in a similar manner that 

oo 

U(F^ l > n ,F?> l > n ) v& llah and lim limsup V U\\F^ l ' n \\l^ = 0. 

l=m 

Next, let us define 

p k>n (j) := carr(A$G, A^. k G), j G Z. 

Assumptions (Al) and (A2) imply (cf. Lemma 1 of [2]) that there exist a sequence 
{PkU)) jeZ such that \p k ,n(j)\ < PkU) for a11 j G Z and n G N, with 

oo 

PkUf < °°- 

j=-oo 

Now, for any s, t G [0,T], we have 

[s/(«„A„)] [t/(«„A n )]-i 
| Z! ^Un )jP 2,/,n^| < na 2„ nAn £ £ KnMi ~ + K/2])|' 

i=[k/u n ]+l j=[k/u n ]+l 

oo 

<C Yl Pk( U nJ + K/2]) 2 -> 
j=-oo 

when n — )• oo, where C" > is a constant that depends on s and t. 

Additionally, for any m = 1, — 1, we estimate the contraction F^' l,n (g> m F^' l ' n via 

II jrhhn <> plA"!^ 

-Tj ||^<g>2(i-m) 

[t/KA„)] 

< (u n A n ) 2 ^ kfc,n(«nKl - «2|)| m kfc,n(«nK3 - M|)| m 

» 1,42,13:14=1 

x |r fc>n (u n |ii - i3|)| i_m |rfc,n(«n|*2 - h\)\ l ~ m , 
the r.h.s. of which is bounded by a positive constant times 

[t/(u n A n )} 

u n A n Y \rkA u nh)\ m \rk,n(u n i2)\ m \rk,n( u n\k - h\)\ l ~ m \rk,n{un\h - h\)\ l ~ m 

[t/(u n A n )} . [t/(u n A n )] . 2 

= u n A„ ^ ( Y \rkA u nh)\ m \r k>n {u n \i 1 -i 2 \)\ l ~ m \ . (A.3) 

12=1 ^ H=l ' 
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By the Cauchy-Schwarz inequality, the r.h.s. of (A. 3) is in turn bounded by 

[t/(u„A n )] . [i/KA„)] 



Y i^M)! 2 " 1 1+ Y i^(vi)i 2(M 

i 2 =l ^ ii=l ' 



[t/(u n A n )} 

<c y 

12=1 



[t/(u„A„)] 



ti=i 



where C" > is a constant that depends on t. Thus, we have \\Ft' l ' n ® m Ff"" ' 2 



as n — > oo. The contraction F. 



2,l,n 



m Ft' l,n can be treated in an analogous manner. 



"^®2(!— m) 



The preceding steps and Theorem 5 of [2] allow us to conclude that the finite dimen- 
sional distributions of the process {Z\, ^t 2 )tg[o,T] converge weakly to those of 



Y Ua f B t = y m ^- m l B t: te[o,T). 

\ 1=2 

Stability of the convergence and tightness in D 2 ([0,T]) can be established similarly as in 
the proof of Theorem 6 of [2] . □ 
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